Cerenkov Emission by Neutral Particles in Gravitoelectro- 
magnetic Fields 
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It is shown that under the post-Newtonian approximation the Einstein 
equations can be reduced to the standard Maxwell-type field equations 
in a medium; in such a context the Cerenkov emission by a neutral 
particle gives large energy loss while the particle moves at faster than 
the phase speed of waves in the medium. 
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It has often been noted that magnetism can be understood as the consequence of electrostatics 
plus Lorentz invariance. Similarly, Newtonian gravity together with Lorentz invariance in a 
consistent way must include a gravitomagnetic field. This is the case of gravitoelectromagnetic 
(GEM) form of the Einstein equations in a medium [1-5]. In such a framework, the refractive 
index n = v^of a medium is often is more than one [see (27), below]; note that the statement 
that only accelerated particles radiate is not correct when there are waves with n > 1 [6] . Hence 
the Cerenkov emission by neutral particles in constant rectilinear motion is unavoidable. 

It is well known that in the slow-motion limit of general relativity, accurate to post- 
Newtonian (PN) order, where the field equations can be reduced as [7]: 
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VV = 47rGp, 



(1) 



V^A = 167rGpv/c; 



(2) 



moreover the harmonic gauge condition and the Lorentz-type force on a unit mass reduced to 

V.A+1|^ = 0, (3) 

dv ^ IdA V ^ ^ 

~_V0--— + - X Vx A, 4 
at c at c 

where (p is the Newtonian gravitational potential, A/c = qiq is the mixed metric and G is the 
Newton's constant. 

One defines the GEM fields via 

E, = -V0--— , VxA = B, (5) 

in direct analogy with electromagnetism; it follows from these definitions the Maxwell-type field 
equations in the medium [8, 9] 

V X B, = l^j, + V-E, = 47rp„ (6) 

1 rf¥K 

VxE, = --^, V-B, = 0, (7) 

where pg — —Gp is " matter density" and jg = —Gpv is " matter current". Then Eq.(4) 
becomes 

^.E, + ^xB,. (8) 

And all terms of 0{v'^/c^) are neglected in the above analysis. Now one decompose the total 
matter density and current in equation (6) : Pg — Pg + Pgo,ig — jg + 3go, in which the Pgo and 
jpo are the parts of a external field source in the medium, corresponding to a local additional 
mass disturbance, for example, by nonlinear interactions between the fields and the medium, or 
some test particles. In this case it is convenient to definite the vector quantity Dg(t, r) through 
( in what follows one neglects the marks ' ~ ' in p^and j^, for simplicity ) 

t 

Bg^Eg + An J dt'ig{t',r)- (9) 



by use of (9) and the mass continuity equation, 



dpcc 
dt 



+ V • = {a^g, go) , 



(6) and (7) are deduced to 



V X B - ' ' ^ 



1 (9B 

VxE^ = — V-B, = 0. 

c at 



(10) 

(11) 



A dimensional analysis of equations (1), (5), (9) and (11) suggests that it is insightful to 
use the cgs units, since these field variables have the following dimensions: 

[D,] ^ [E,] ~ [B,] ^ [V0] ^ [v]/s, 

thus, one has 



El/G ~ B'/G ~ pE'/pG 



pv 



and by substituting 



E = -Eg/AVc, B = -Bg/4VG, D = -D^/4\/G, 
j = -ig/VG^y/Gpv, p^-pJy/G^VGp, 



(12) 
(13) 



the GEM field equations can be expressed in the standard form 



„ „ 47r. 1(9D „ ^ 19B 
VxB = — J0 + - — , VxE = — — 
c c at c at 



V • D = npo, V • B = 0; 



(14) 
(15) 



then Eqs.(9) and (8) become 



D = E + TT J dt'i{t',r) 



(16) 



o 



at 



E+ - X B 

c 



(17) 



To consider the responses of the medium on the GEM fields, one must introduce material 
relation, which describes the GEM properties of the medium. In view of (16), the states of the 
medium not only depend on a given time-space point {t, r) but also depend on previous times 
and at any point of the medium. Hence, by general reasoning in physics ( independent of a 
specific model for the medium ) one can state that this is a no-local linear relation in the limit 
of linear response ; whose Fourier representation is 

ji{uj, k) = aij{uj, ^)Ej{uj, k). (18) 

In special, for common continuous medium, where the " spatial dispersion " (dependency on 
k) is not important , the relation (18) is deduced to (uy — > crSij) 

j(a;) = a(a;)E(a;), 

i.e. an " Ohm's gravitational law " [10] , where a is "gravitational conductivity" [11] ( for 
details see Refs [10,12]) . Then by use of Eqs.(16) and (18), one has 

Biioj, k) = Eijioj, k)Ej{u;, k), (19) 

here Sij, called the "dielectric tensor", is determined as 

eij{uj, k) = 6ij H crij{uj, k) (a; 7^ 0). (20) 

In fact , if we consider the kinetic description for a self-gravitating system, on basis of Vlasor 
equation of the gravity [13,14], in which the force is Lorentz-type one [Eq.(17)], one can expand 
the distribution function / in powers of the field E; and by Fourier transformation 

A(r, t)= J A{uj, k)e-^('^*-^-"-)da;dk, A{uj, k) ^ J A(r, t)e'^'-'-''-''^ (21) 



A 



it yields the material relation (18) for the first order current j oc / v/i(r, v, i)o?p , in which 
aij{u!, k) is relative to the particle distribution in the medium . 

By applying the Fourier expansion, the field equation (14) becomes 
z(k X B)i = — -£ijEj + —jo,i 

UJ 

kxE = -B; 

c 

by elimination of B, 

47ri 

Aij{u;, k)E^{u, k) = jo,i{u;, k), (22) 



where 

A,,-(a;, k) = ^(^ - 6,j) + e,j{u;, k). (23) 
Multiplying Eq.(22) by any modular vector e^*^, with €^ ^6^*^ = 1, one obtains 

(e-^el\El'^^^u;[ei*.f{u;M , (24) 



in which 



4 ^ £:,k = ^U^: ^K,iei:j + ^(k • e^)(k • eD; (25) 



UJ 



and £^is the "dielectric susceptibility". By neglecting the terms in right hand ( the nonliear 
terms), one gets from Eq.(24) the homogeneous wave equation; and from here one finds the 
dispersion equation 

1,2 2 

el-'-^-O . (26) 



Let UJ = cc;'^(k) be its a solution ; physically, the progressive wave in the opposite direction is 
also possible, thus Eq.(26) has another solution with negative frequency uj — — a;'^(— k). 



Consider the lowest order equation of motion for a test particle 



for plane waves oc e «'^*+«kr^ one has : — ia;Vk,a) = — 4-\/GEk,a;, thus 

fCo 



—icu 



I.e. 

G p 

jk,i = CTijEkJ, CTij = CT^ij- = 4— 

For transverse mode (k • Ck = 0), on the other hand , one has from Eq.(20) and Eq.(25) 

4 = l + ,_ = l + _/; (27) 

with the "gravitational conductivity" a, 

Hence from Eq.(26) we obtain the dispersion relation 

o> F> ^28) 

uj ^ kc {cu cUq). 

To calculate this emission radiated by the moving external currents, we form the scalar product 
the first equation of Eq.(14) with E, of the second one with B and take the difference; then 
integrating over a volume and using the Gauss theorem, as a result we obtain 



taking into account that the E and B vanish at infinity, thus 

= _ E- — + B • — = -jo • E. (29) 
dt 4.1:^ dt dt^ ^ ' 

i.e., the work done on the gravitoelectric fields per unit time, by the external currents, makes 
the GEM field energy density ( [jqE] ~ [jogEg/ G] ~ [pv^] / s) get increments. Hence the GEM 
wave energy radiated is 

U'' ^- J dtj dTE{r,t)-ioir,t). (30) 

Using the Fourier transforms of E and jo and taking account of the condition of real field, we 
obtain the expression 

U" ^- (27r)^ Re J dudki^ {u, k) • E* {u, k) . (31) 

In view of Eq.(24), 

U' = 2 (2.f Re I rf^rfi,'"'^-^|e"(k)^J. (^.k)r 

by use of the Landau rule and the Plemelj formula [15], 

Im-, — = Im 2 = —T^^ i T^l) 

one has 

= {'^T^fj P dwujc-^\e''*{}^) •j(a;,k)|^(5 (^"^ - '^el^ 



where 



s^ = <k = 4(^'k)e^,,e^*; (33) 



the second term of Eq.(25) has no contribution to the integration in Eq.(32); and in the last 
step we consider the solutions lo — ±a;'^(k) of the dispersion equation, A;^ — = 0, have the 
same contribution to the integration. 

On the one hand , one may consider any instantaneous velocity as uniform motion at 
short time intervals , and on the other hand the space should be approximately fiat for the slow- 
motion and weak gravity, hence it is reasonable to study the emission radiated by a neutron 
particle with mass m in constant rectilinear motion passing through the medium. Then the 
current density is 



where go = \fG'm, r{t) is the position vector of the motion particle. Considering that r = 
ro + vt, V = const., the Fourier transform of Eq.(34) is 



Because of the delta-function in Eq.(35), for the emission by a uniformly moving neutron particle 
has to be satisfied the Cerenkov condition: u;"' — li-v, or v > u;"'/k — c/n{uj), which essentially 
corresponds to the laws of energy and momentum conservation in the radiation processes. In 
view of Eq.(27), then n{uj) > 1; therefore, the Cerenkov emission by the neutral particles in 
constant rectilinear motion is unavoidable. 

Substituting the equation Eq.(35) into Eq.(32) the foregoing results yield readily 



jo(r,t) =gov(t)5(r-r(t)). 



(34) 



jo(t^,k) 



(27r)3 



vexp(ik • ro)S{u! — k • v). 



(35) 




(36) 



Using the following identity of delta-function 




and considering the definition of radiated power, I"^ — lim —we obtain 




(37) 



o 



For transverse waves the polarization vector e*(k) in Eq.(37) has two independent components 
(eW,e(2)), then the total power radiated is obtain by adding the two parts, i.e., proportional to 



.(1) 



(k)-v 



+ 



e(^)(k) ■ V ) ; and because e*^^\ e*^^^ and k(= k/Zc) constitute orthogonal frame 

Sij — kikj/k'^, one has 



in the complex space, e\^^ef'^* + ef'^ef^* 



5«(k) 



+ 



3(2) (k) 



v^-(k-v)^ = 'i;^sin2^; 



hence 

= J^Jl-^\ (38) 

It is well known that the gravitational system can emit energy in quadrupole radiation; and 
the power radiated is proportional to the order of c~^. But one see from Eq.(38) that the 
Cerenkov emission occurs in the order of c~^. Therefore, in general, the Cerenkov emission 
is more than the quadrupole radiation. For example, consider a binary system, in which 
the accreting compact object is a larger black hole with mass M = 10^ Mq and companion 
star with mass m = Mq moves in nearly circular orbit around the black hole at the radius 
oo = 300rg {vg — 2GM/c^, the Schwarzschild radius ) . In such a case, from the Kepler 
law, (27r/p)^ = G{M + m)/al, then one gets p = 1.56 x 10^(/i); at the same time, we may 
taket; ~ (2^7^/00)^/^ ~ 1.6 x 10^{ cm ■ s ). Furthermore, one can estimate the medium 
density at r = ao, p = 2.9 x 10~^(ao/3rg)^/^ = 2.9 x 10~^{g ■ cm~^) [16]. As a result, due to the 
quadrupole radiation, the energy loss for the system is [17] 

— = 3.0 X 1033(-^)2(il)4/3(^^)-10/3 ^ ^ ^ ^ IQ^^ ( . ^-1) 

here n = Mm/ (M + m). And one gets from Eq.(38) the energy loss in Cerenkov radiation 
7' = 3 X 10^^ { 7) i ~ 4.8 x lQ^\erg ■ s'^)- 

we have taken co 10~'''(s~^) at the above. Then = £* ~ = 2.5 x 10^, thus the Cerenkov 
condition is satisfied: v > c/n{uj). 

n 
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